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^_^ Using minimal models for low Reynolds number passive and active rotors in a fluid, we characterize the hydrodynamic interac- 

^-H tions among rotors and the resulting dynamics of a pair of interacting rotors. This allows us to treat in a common framework 

^^ passive or externally driven rotors, such as magnetic colloids driven by a rotating magnetic field, and active or internally driven 

vN rotors, such as sperm cells confined at boundaries. The hydrodynamic interaction of passive rotors is known to contain an az- 

^ imuthal component ~ 1/r^ to dipolar order that can yield the recently discovered "cooperative self-propulsion" of a pair of rotors 

/-^ of opposite vorticity. While this interaction is identically zero for active rotors as a consequence of torque balance, we show 

^^ that a '^ 1/r^ azimuthal component of the interaction arises in active systems to octupolar order. Cooperative self-propulsion, 

^O although weaker, can therefore also occur for pairs of active rotors. 
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There has been a lot of interest in the properties of individ- 
ually propelled particles that take up energy from their en- 
vironment and collectively generate motion and mechanical 
stresses^. These systems exhibit remarkable collective behav- 
ior, quite different from their equilibrium counterparts. One 
important class of propelled particles that has received rela- 
tively little attention is "rotors", that is particles that rotate in 
a fluid in response to externally applied or internally generated 
I torques. An external torque is provided for instance by a rotat- 
ing magnetic field ^^ applied to particles carrying a magnetic 
moment, or by an optical trap^ ^. Internally generated torques, 
on the other hand, originate from forces exerted by the rotor it- 
self on the fluid. We will refer to the first class of rotors driven 
by external torques as passive rotors, while the name of active 
rotors will be reserved to those that are internally driven. Sev- 
eral examples of active rotors are found in the living world, in- 
cluding sperm cells'^^^, bacteriat^SC^and algae^near a solid 
surface. Various artificial swimmers, inspired by their living 
counterparts, have also been engineered over the past decade, 
and provide realizations of active rotors ^^"^^^ . Other examples 
of internally driven or active rotors include the rotating motors 
found at the basal bodies of cilia and flagella anchored at the 
cell membrane— and ATPase molecular motors embedded in 
fluid membranes'^. In some cases, it is also possible to ap- 
ply an additional external drive to an active swimmer, as with 
magnetotactic bacteria in a rotating magnetic field ^^. 
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Rotors generally move through a viscous fluid and gener- 
ate flow in the fluid. This flow in turn affects the dynamics 
of other nearby rotors. The goal of this paper is to charac- 
terize the flow field induced by various types of rotors and 
the resulting hydrodynamic interactions among them, high- 
lighting the differences between active and passive rotors. We 
then discuss how the flow field controls the rotors' dynamics. 
We limit our analysis to the case of rotors confined to a two- 
dimensional plane with their spinning direction perpendicular 
to that plane. This constraint is often realized in experiments 
by confining rotors to a liquid/air interface^ ^, or a thin layer 
of viscous fluid embedded in a different unbounded viscous 
fluid ^^ ^^, or in the case of many swimming unicellular or- 
ganisms simply because they are attra cted to a boundary that 
in turn triggers the rotational dynamics '^'^ ' ^^"^^ '. 

The main difference between passive and active rotors re- 
sides in the nature of the azimuthal component of the flow 
they create. Passive rotors transfer external torque to the fluid, 
generating a dipolar azimuthal flow that decays as 1/r^, with 
r the distance from the rotor. Active rotors do not transfer any 
net torque to the fluid, as the torque transmitted from the rotors 
to the fluid is balanced by an equal and opposite torque trans- 
mitted by the fluid to the rotors. As a result, the 1/r^ dipo- 
lar azimuthal component of the flow field is zero for active 
rotors. This important distinction has been known for some 
time and was first pointed out in a classic paper by Batche- 
lor^^. More recently, it has been stressed by Lenz and collab- 
orators ^^^^ and by Michelin and collaborators^^, but is still 
sometimes overlooked in the literature. Our work shows that 
while the 1/r^ azimuthal component to the flow field vanishes 
as expected for active rotors, the cycle-averaged flow field of 
these force-free swimmers does contain a non-zero azimuthal 
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in-plane component, although only to octupole order in the 
multipole expansion. The corresponding contribution to the 
flow field hence decays as 1/r^ rather than as 1/r^ as for pas- 
sive rotors. 

This finding is an important new result of our work be- 
cause the existence of a finite azimuthal interaction alters sig- 
nificantly the collective dynamics of the rotors. It is known 
that an azimuthal interaction between two passive rotors of 
the same vorticity (i.e., spinning in the same direction) re- 
sults in the pair revolving around each other along a circular 
path, as seen in experiments on externally driven rotors^ ^. It 
was shown recently by Leoni and Liverpoolf^^that for two ro- 
tors of opposite vorticity the same 1/r^ azimuthal interaction 
turns the rotor pair into a self-propelled unit, in that the two 
rotors push each other in a common direction. The assump- 
tion of a 1/r^ azimuthal interaction made in Ref. 27 seems 
to restrict the result to the case of externally driven rotors, al- 
though much of the discussion in that paper refers to active ro- 
tor. Our work demonstrates that this type of "cooperative self- 
propulsion" also occurs for active rotors and arises not from a 
~ 1/r^ azimuthal coupling, that vanishes for force-free parti- 
cles, but from the new 1/r^ octupolar contribution to the flow 
that we have identified. We also show that for active rotors the 
additional attractive and repulsive interactions arising among 
pullers and pushers from the radial dipolar component of the 
flow field can lead to the formation of a bound self-propelled 
pair consisting of two oppositely rotating swimmers. 

To highlight the difference between passive and active ro- 
tors, in the following we first review known results for the sim- 
plest passive rotor consisting of a sphere driven by an external 
torque. Then we analyze the flow field and hydrodynamic in- 
teractions of a minimal model for an active rotor, consisting a 
one-bead swimmer that rotates by pushing on the fluid with a 
force oriented at a finite angle to the swimmer's long axis (see 
Fig.|2j. 



given by 



^ii^) = ^Jd^'^iji^-^')M^')^ 



(2) 



where ^ij{r) is the Green's function of Eqs. ([T]) and depends 
on the geometry of the problem. In the following we use the 
Oseen tensor given by ^ij{r) = ^{Sfj -\-rifj) with ff = ri/r, 
which corresponds to an unbounded fluid in three-dimensions, 
but can also describe the flow field generated by particles 
confined to a fluid-air interfaceE^EII in the plane of the in- 
terface. Our results are, however, easily extended to other 
Green's function describing different physical situations (see 
section [4]). 

The rotors dynamics is controlled by force and torque bal- 
ance, given by 



F^+F^ = 0, 







(3a) 
(3b) 



where the superscripts d, a and e stand respectively for drag, 
active and external, and the forces and torques are exerted on 
the rotors. We only consider below situations where the exter- 
nal force is zero (in particular, the rotors are assumed to be at 
neutral buoyancy), but generally allow for the presence of an 
external torque. 



2.1 Single passive rotor 

We model a passive rotor as a single neutrally buoyant sphere 
of radius a rotating under the action of an externally imposed 
torque, T^ = T^z. In this case there are no active forces or 
torques, and the drag force and torque exerted by the fluid on 
the sphere are given by 



2 Model 

We consider rotors moving through a viscous incompressible 
fluid in the xy plane and spinning around the z axis. Due to 
the small size of the rotors, inertial effects are negligible both 
on the dynamics of the fluid and that of the rotors. The fluid is 
therefore described in the zero Reynolds number limit where 
the flow field u is solution of the incompressible Stokes equa- 
tion 



T]V^u-V/? = -f(r) 
Vu = 0, 



(la) 
(lb) 



where r] is the viscosity of the fluid and f is the force density 
exerted by the rotors on the fluid. The solution of Eqs. ([T]) is 




(4a) 
(4b) 



where ^ = 67tr\a and ^^ = ^7tr]a^ are the translational and 
rotational friction of the sphere, re is the position of its cen- 
ter, Vc = dvc/dt is its velocity and fi) is its angular velocity 
describing rotations about the center. The flow u^(r) is the 
one that would be obtained if the sphere was removed and re- 
placed with fluid. The laplacian term on the right hand side of 
Eq. ( [4a| is the Faxen law correction for the friction on spheri- 
cal objects in a viscous fluid and is important when flow gra- 
dients occur on length scales comparable to the sphere's diam- 
eter. ^^ The translational and rotational velocity of the sphere 
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are immediately obtained as 



Vc 



i + ^vMu^ 



J re 



(0= -i^z+- [Vxu'^jj 



^R 



(5a) 
(5b) 



If there is no externally imposed flow and only one rotor is 
present in the fluid, then u^(r) = 0. In this case a single pas- 
sive rotor does not translate (vc = 0) and it rotates around the 
z axis with angular velocity G) = zT^/^^. 

The flow field due to the rotation of the sphere is obtained 
by enforcing a no- slip boundary condition at its rotating sur- 
face and is given by 122 



u(r) 



Snrfr^ 



rz X r . 



(6) 



The flow field is purely azimuthal in the xy plane and it decays 
as 1/r^ (see Fig. T|. 



where ^^^ = (5^^^ - 3rjrk)ri/r^ and Wtjk = £mq£njkrq/r^ are 
the jk symmetric and antisymmetric parts of dj^^ij. The 
decomposition of the flow into its Jifijk part (known as a 
stresslet) and ^tj^ part (known as a rotlet) coincides with a 
polar decomposition in the plane perpendicular to the axis of 
rotation: the stresslet is purely radial while the rotlet is purely 
azimuthal.The antisymmetric part of the integral represents 
the total torque exerted by the rotor on the fluid. Torque bal- 
ance requires that this equals the external torque applied to the 
rotor, or 

eikjjdr'r'jj{r') = xt. (8) 

As a result, the antisymmetric or rotlet part of the flow field 
given in Eq. ^ equals the azimuthal flow field given in Eq. ^ 
for all passive rotors and vanishes when T^ = 0. In other 
words, the azimuthal flow field of an externally driven rotor 
is to leading (dipolar) order identical to that of a sphere driven 
by the same torque, regardless of the details of the rotor. The 
symmetric part of the integral, on the other hand, is not con- 
strained and depends on the geometry of the rotor. 
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2.2 Single active rotor 

An active swimmer is one that propels itself by pushing on 
the fluid. In order to obtain net propulsion in the absence of 
an externally applied force or torque, the drag center C and 
the thrust center T of the swimmer must be located at two 
different points. With this in mind, we consider the simplest 
realization of a rotating swimmer, consisting of a spherical 
bead of radius a centered at point C and providing the drag and 
a point force — F exerted on the fluid at the end T of a fictitious 
"flagellum" of length i providing the thrust. We denote by V 
a unit vector pointing along the swimmer's long axis (chosen 
as the direction of the flagellum), and by the angle between 
V and F (see Fig. [2]). The dynamics of the swimmer in the 



Fig. 1 (color online) Map of the flow field generated in the xy plane 
by a passive rotor modeled as a sphere driven by an external torque 
(represented by a circular arrow inside the rotor) directed along the 
out-of-plane direction z. The rotation is counter-clockwise and 
occurs around the z axis. The fluid velocity is shown as arrows with 
color coded magnitude. 

It is important to note that this azimuthal flow field is a gen- 
eral property of passive rotors. This can be seen by expand- 
ing Eq. ([2]) in multipoles. The monopole term vanishes be- 
cause there are no external forces acting on the rotor. The first 
nonzero term in the multipole expansion is the dipolar term 
that has the form 

^Kr) = -^l^ijkir) ^<^ijk{r)]Jdr^r[fj{r') , (7) 



V" 



-F 



^C 




Fig. 2 Left: Schematic of the single-bead swimmer. The swimmer 
exerts an active force — F on the fluid applied at point T at an angle 
to the swimmer's axis defined by the unit vector V. The bead has 
radius a and its center C is the hydrodynamic center of the swimmer. 
The rotation is counterclockwise about the z axis normal to the page 
through C. Right: Schematic of the swimmer's motion. The 
swimmer, drawn at three different times, undergoes a full rotation 
around itself every Tq = In/coo while translating with velocity vq, 
resulting in a circular trajectory of radius Rq = \()/(0(). 
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absence of external torque is governed by Eqs. ^ with T^ = 
and F^ = F, while the active torque T^ = —iv x F is equal to 
the moment in C of F^. The drag force and torque are given by 
Eqs. (|4]). For an isolated swimmer, u^ reduces to the flow field 
created by the thrust force — F exerted on the fluid at point 
T. The linear and angular velocities of the swimmer are then 
given by 



Va 



3a 



a-" 

4^3 



{dij^ViVj)---^{5ij-3ViVj) 



(0 = 






VxF 1 
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(9a) 



(9b) 



The thrust force — F exerted by the swimmer on the fluid at T 
affects the dynamics in two ways: (i) "directly" through the 
reaction force +F exerted back on the swimmer by the fluid 
at r, and (ii) "indirectly" through u^. These two contributions 
cancel each other when £ = a, i.e., when the thrust force is 
exerted at the surface of the bead, giving Vc = and fi) = 
in this limit. Conversely, the "direct" contribution is dominant 
when a<^L Since V and F both lie in the xy plane, fi) = (Oqz 
is along z- An isolated swimmer then translates at Vc while 
rotating at an angular frequency (Oq about an axis through its 
center C directed along the z direction. The trajectory of C is a 
circle of radius R^ = vo/ftio in the xy plane, as shown on Fig.[2J 
where vq = |vc|. The period of rotation is Tq = 271 /cOq. Note 
that when (j) ^ 0, Rq diverges and the trajectory becomes a 
straight line, as expected for a linear swimmer exerting a trust 
force directed along its long axis. 

In order to keep the model tractable, we work in the a <C ^ 
limit. Since Rq/I ^ ^ (^{a/i)^) /sin0, the radius of the circle 
is then negligible for all finite values of (j) and the swimmer is 
merely rotating rigidly about its center K Moreover, the dy- 
namics is dominated by the "direct" contribution of the thrust 
force and the velocity and angular velocity are simply given 
by \c^^/l^ and (D c^ -^Fsin^z/C^. 

The flow field created by the one-bead swimmer is the sum 
of the flows produced by three sources: a point force — F lo- 
cated at T corresponding to thrust, the translation of the bead 
with velocity v = F/^, and the rotation of the same bead with 
angular velocity G) = —£v x F/^^. This gives 



Ui{r) 



^ij{r-rc^iv)Fj^- ^ij{r-rc)Fj 



Sttt] 



Sttt] 



where we have used Ft = re — iV and expressed v and (0 in 
terms of £, V and F. In order to study the the far field behavior 

t A detailed analysis of the role of a finite a, including the consequence of Rq 
not being neghgible, can be found in the Electronic Supplementary Informa- 
tion. 



of the flow, we expand Eq. ([TO]) in powers of £/r. The flow 
field can then be written in the form of a multipole expansion. 



as 



u(r) 



I" 

n=0 



(n) 



(r), 



(11) 



where u*^")(r) ~ [F /r[r){l/rY represents the (n + l)-pole 
contribution to the flow. In the absence of external forces the 
monopole term u(^''(r) vanishes. The lowest non- vanishing 
term of this expansion is the dipolar one. We note, however, 
that as a consequence of torque balance the rotlet due to the 
rotation of the bead, given by the last term on the right hand 
side of Eq. ([TO]), is precisely canceled by the rotlet component 
of the first two terms, representing the flow due to the force 
dipole formed by the two point forces. As a result, only the 
stresslet (radial) part of the flow is nonzero to dipolar order, 
with 



«P(r) = 



1 



IF 



iFjVkJ^ijk 
cos0 + 3cos(20- 



(12) 



where 6 is the angle between v and r. As we pointed out in 
section |2.1[ this is a general consequence of torque balance 



and holds for any rotating swimmer. 

Retaining terms up to third (octupolar) order in the expan- 
sion, the total flow field of a single one-bead swimmer is given 
by 
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Sttt] 



^Vk^i, 



ijk ' 



-h^idkl^ij 



-^VkVl^mdklm^ij 



Fj^e 



, (13) 



where all the derivatives of the Oseen tensor are evaluated at 
(r-rc). 

On time scales longer that the period of rotation Tq of 
the swimmer(see Fig. [2]), it is useful to consider the cycle- 
averaged flow field, denoted by ( ...) below. As shown in Ap- 
pendix |A| the average over a period can be obtained by aver- 
aging the instantaneous flow over the orientation v(^) of the 
swimmer's axis. Since the active force Fj in Eq. ([13]) can be 
written as Fj = Fcos v^ — F sin0 (z x v) •, it is then evident 
that the dipolar and octupolar terms (first and third term on 
the left hand side of Eq. ([13])) contain even powers of the unit 
vector V, while the quadrupolar term is cubic in v. The lat- 
ter will therefore vanish when we average over the direction 
of the swimmers' axis. This result is generic for swimmers 
confined to a plane. It should be stressed, however, that the 
instantaneous quadrupolar flow field is of course nonzero at 
all times. 
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and neglecting ^{a/tj corrections to all the coefficients, the 
total rotation-averaged flow of a single one-bead swimmer is 
given by 



(u(r)) 



Fcosij) 



^2 



3Fsin0 t 



COS0 f 



— rzx r- 



^ . (16) 



The azimuthal term is proportional to the active torque ^F x 

V = — ^Fsin^z, while the radial term is proportional to F • 

V = Fcosij). The latter is zero for a "pure rotor" (F • V = 0), 
attractive for a "puller" (F • V < 0) and repulsive for a "pusher" 
(F-V>0). 



Fig. 3 Map of the instanteneous flow field generated in the xy plane 
by a one-bead pure active rotor (i.e. thrust force perpendicular to the 
main axis) rotating counter-clockwise around the out-of-plane 
direction z. The fluid velocity is shown as arrows with color coded 
magnitude. The arrows represent the forces (straight arrows) and 
torque (curved arrow) exerted by the rotor on the fluid. The inset 
shows the flow field averaged over a rotation. The averaged flow 
field is azimuthal and dominated by the average flow field of a point 
octupole. Distance is normalized by the length i of the rotor. 



The cycle-averaged dipolar stresslet contribution to the flow 
is isotropic and radial and is entirely controlled by the compo- 
nent F cos of the propulsive force parallel to the axis of the 
dumbbell, with 



i(i) 



(r) 



IF cos ^ 



(14) 



The quadrupolar term in Eq. ([13]) changes sign when the ro- 
tor is rotated by 71 and thus does not contribute to the cycle- 
averaged flow. The octupolar term gives rise to both a radial 
and an azimuthal contributions upon cycle averaging, with 



,(3) 



(r) 



i^F 



167rr]r4 ^8 



- cos 0f+-sin0zxf 



(15) 



The radial part will be neglected below as it yields a correction 
to the cycle-averaged dipolar flow, with a faster decay. The 
azimuthal part, on the other hand, represents the leading order 
azimuthal flow. Both Faxen effects and the finite value of the 
radius Rq of the circular trajectory of the center of the rotor 
give corrections of order {a/i)^ or higher to the coefficients of 
both the radial and azimuthal octupolar flow fields, but do not 
change their qualitative behavior and will be neglected here. 
The full expressions for the cycle averaged octupolar flow field 
can be found in the Electronic Supplementary Material^ . 

Keeping only the leading contribution of each component 
(i.e. the radial dipolar term and the azimuthal octupolar term). 



3 Hydrodynamic interaction of two rotors 

In the presence of other rotors, the flow field generated by one 
rotor yields hydrodynamic interactions on the others. These 
interactions in turn control the collective rotor dynamics. To 
understand the role of such hydrodynamic couplings, in this 
section we consider the effect of the fluid mediated interaction 
on the dynamics of pairs of rotors. Again, we will use mini- 
mal models for both passive and active rotors, corresponding, 
respectively, to a single rotating sphere and to the one-bead 
swimmer discussed in the previous section. 

We consider two rotors in the xy plane. The positions of 
their drag centers are denoted by vectors {Ya}a=i2- ^^ addi- 
tion to the fluid mediated hydrodynamic interactions, the two 
rotors are assumed to experience a short-range central repul- 
sive force that opposes their overlapping. The equations of 
motion of the two such rotors have the form 



dtra 



ra 



W 



I5)^a(5 



(17) 



where up is the flow field induced by the rotor j3 and g{ro^p) 
is the magnitude of the repulsive central force and vanishes 
at distances larger than the effective diameter d of the rotors, 
with r^^^ = r^ - r^ ^ . 

For passive spheres, the flow field u that appears in Eq. ([17]) 
is given by Eq. ([6]). For the active one-bead swimmer, we only 
consider the dynamics on time scales large compared to the 
period of rotation Tq. As shown in Appendix [A| in this limit 
we can replace the instantaneous flow field u in Eq. (Vl) by its 
cycle averaged value given by Eq. ([16]) ^ . This can be written 

J When an explicit form was needed for the repulsion force (e.g., to run sim- 
ulations, not shown), we used g{r) = [{r/d)^ — l] /r^ if r < d and g{r) = if 
r> d. 



t If the radius Rq ^ 



I of the trajectory of a single rotor were not negligible. 



the fact that both rotors are moving around a circle would have to be taken into 
account when performing the period average. This case is treated in detail in 
section 2 of the Electronic Supplementary Information. 
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in a general form that applies to both passive (representing in 
this case the instantaneous flow) and active rotors as 

(u^(r«)) =F«'^(r«^)r«^+A«'^(r^^)e^zxr^^ (18) 

The first term is the radial flow, with F^ = 
(^Fcos0)/(16;r77r^) and F^ = 0. The second term is 
the azimuthal flow of strength A, with A^ = T^/(8;rr]r^) and 
A^ = {3fFsm(j))/{64nrir'^). A parameter £a = ±1 is intro- 
duced to characterize the vorticity of the rotor, with £a = ^l 
corresponding to clockwise (+1) and counter-clockwise (-1) 
rotations. 

The term proportional to a^ in Eq. ([17]) contributes to the in- 
teraction to octupolar order through the laplacian of the dipo- 
lar flow field. For a passive sphere, the dipolar flow is a rotlet 
and its laplacian is zero. For a one-bead swimmer, the dipolar 
flow is given by Eq.p^and its laplacian has the same structure 
as the radial part of u^\ with an extra (a/i)^ factor. We can 
therefore neglect it as well. 

To study the coupled dynamics of the rotor pair we in- 
troduce the center of "mass" and relative coordinates, R = 
(ri + r2)/2 and r = ri — r2, respectively. Without loss of gen- 
erality we assume that rotor 1 is rotating clockwise £i = +1. 
The equations of motion of the two rotors can be rewritten in 
the form 



(9,R=-A^'^(r)(l-e)zxf , 



(9,r = A^'^(r)(l + e)zxf + 2 



F^^P{r) 



Kr) 



(19) 
r , (20) 



where e = +1 corresponds to "like" rotors, rotating in the 
same direction, while £ = —I describes "unlike" rotors, with 



opposite circulations. For rotors of equal vorticity Eq. (19 ) for 
the center of mass coordinate reduces to dtR = 0, indicating 
that in this case the position of the center of mass of the ro- 
tor pair remains fixed at its initial value. We now discuss the 
dynamics of various pairs of active and passive rotors. 



(a) 



(b) 



o 





Fig. 4 Schematic trajectories of two identical rotors of equal 
vorticity. The two dots are the initial positions of the rotors. A line 
with alternating color is used where the trajectories are 
superimposed, (a) Pair of passive rotor or pure active rotor (circle), 
(b) Pair of contractile swimmers (shrinking spiral), (c) Pair of 
tensile swimmers (expanding spiral). 



1. Like passive rotors. In this case £ = 1 and the center 
of mass is static. The radial component of the force re- 
duces to the short-range repulsion force. The two rotors 
rotate about their center of mass with an angular velocity 
Q. = T^/{47irirQ) describing a circle of diameter tq. The 
diameter ro is simply the initial separation of the rotors if 
they start out at a separation greater than the range d of 
the repulsive interaction and equals d if the initial sepa- 
ration of the centers of the two rotors is less than d. The 
trajectory is as shown in Fig. (|4^). 

2. Like active rotors. Again £ = I and the center of mass of 
the pair remains fixed. The rotors experience active ra- 
dial forces F^ that can be attractive (0 < \(j)\ < ;r/2, cor- 
responding to pullers) or repulsive (7r/2 < |0 1 < tt, corre- 
sponding to pushers). The relative coordinate of the pair 
of swimmers describes a circular trajectory of varying ra- 
dius, corresponding to an inward bounded spiral (Fig.|4]3) 
or an outward unbounded spiral (Fig. [4]:), depending on 
the nature of the swimmers. Writing r in polar coordi- 
nates as r = (r, 0) we can solve for the dynamics of the 
angular coordinate, which is given by 



0{t) = -ftSin(l) 



To / 3 I 3Ficos^ \ 



For tensile swimmers (pushers) the angular dynamics 
slows down as t~^^^ in the asymptotic regime, while 
for contractile swimmers (pullers), the dynamics is reg- 
ularized by the onset of the short-range repulsion. Then 
6{t) ^ t Sit long times and the two rotors eventually set- 
tle into circling each other describing a circle of diameter 
re where re < d is such that the attractive force is bal- 
anced by the short-range repulsive one. In the limit case 
when (j) = 71 /2 (corresponding to pure rotors), the radial 
interaction reduces to the short-range repulsion and the 
dynamics is similar to that of a pair of like passive rotors 
with an angular velocity Q now equal to A'^/(2r) ^ l/r^ 
instead of A^/(2r)-l/r^ 

3. Unlike passive rotors. In this case £ = —1 and the dy- 
namics of the relative coordinate becomes purely radial 
and, because the rotors are passive, repulsive and short- 
ranged. The center of mass coordinate acquires a veloc- 
ity in the direction orthogonal to the line joining the two 
rotors. The pair of rotors becomes "self-propelled", in 
that it develops sustained motion in the direction orthog- 
onal to the line joining their centers. For initial separa- 
tion ro > d, the two rotors move along parallel lines at 
a speed V = T^/(8;rr7rQ), as shown in Fig. pa. For ini- 
tial separations ro < d there is an initial transient associ- 
ated with the two rotors moving apart due to the repulsive 
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force, after which they settle into the self-propelled tra- 
jectory along parallel lines shown in Fig.|5^ with a speed 

V = t'/{S7irid^). 

4. Unlike active rotors. As in the case of like rotors, the 
trajectory of unlike active rotors is similar to that of pas- 
sive ones, with additional dynamics that arises from the 
repulsive/attractive radial part of the hydrodynamic inter- 
action. For the case of attractive interactions, i.e., pullers, 
there is an initial transient during which the two rotors 
approach each other until the attractive interaction is bal- 
anced by the short-range repulsion. After this time, the 
two rotors move as a self-propelled pair along parallel 
straight lines with a velocity that scales as the inverse 
fourth power of their separation (see Fig. [5})). To un- 
derstand the behavior of pushers it is useful to explicitly 
integrate the equation of motion to obtain the time evo- 
lution of the magnitude of the center of mass coordinate, 
given by 



R{t) = -£^tan0 



^0 



('•o- 






It is then evident that the "self-propulsion" speed of the 
pair decays with time. The active force driving the self- 
propulsion scales as 1/r^ while the force driving the re- 
pulsive separation between the rotors scales as ^ . As a 
result, the repulsive interaction eventually wins and the 
trajectories of the two rotors diverge from each other, as 
shown in Fig. [5]:. As for like rotors, pure active unlike 
rotors (0 = 7t/2) behave like passive unlike rotors, but 
move with a speed V ~ 1/r^ instead of 1/r^. 





(a) 






(b) 
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Fig. 5 Schematic trajectories of two identical rotors of opposite 
vorticity. The two dots are the initial positions of the rotors, (a) Pair 
of passive rotors or pure active rotors (straight parallel lines), (b) 
Pair of contractile swimmers (converging lines), (c) Pair of tensile 
swimmers (diverging lines). 

Summarizing, like rotors circle around each other, while 
their center of mass remains stationary. Unlike rotors are pro- 
pelled in the direction orthogonal to the line joining their cen- 
ters. The difference between passive rotors and active rotors 



primarily arises from the presence of active radial hydrody- 
namic interactions that can be attractive or repulsive depend- 
ing on the nature of the swimming mechanism. When the ac- 
tive propulsion mechanism results in pure rotation (i.e., when 
(j) = 7r/2), the difference between active and passive rotor tra- 
jectories lies only in the scaling of the velocities as a function 
of rotor separation. 

4 Discussion 

The important role of the azimuthal hydrodynamic interac- 
tions between colloidal particles in a fluid driven to rotate 
by an external drive in controlling the dynamical assembly of 
such artific ial microswimmers has been recognized for some 
timel2n3SI, Yet, only recently it has been pointed out that this 
interaction yields a net finite velocity for the center of mass of 
two rotors of opposite vorticity, which will therefore behave as 
a self-propelled pair^^. Such a "cooperative self-propulsion" 
is a novel effect with potential applications. Experimental re- 
alizations of passive rotors naturally generate, however, rotors 
of the same vorticity. The most common realization of collec- 
tions of microrotors consist of magnetic dipoles driven by an 
external rotating magnetic field^^. The direction of rotation 
of the rotors is then imposed by that of the magnetic field, and 
is the same for all the rotors. In Ref.|5l the drive is applied in 
a more local way by using a small coil that is positioned right 
above the rotor. Using two such coils would allow to apply 
opposite drives on two rotors. Similarly, a laser trap can be 
used to control the rotation of an individual rotor ^ ^. Still, ob- 
serving self propulsion with either of these two devices would 
require that the sources of the angular drive follow the rotors 
in their motion. 

Internally driven rotors, on the other hand, do not suffer 
from this limitation and may be a good candidate for the obser- 
vation of cooperative self propulsion. In these active systems, 
however, the internally generated torque that is transmitted by 
the swimmers to the fluid is exactly balanced by the torques 
transmitted by the fluid to the swimmers. As a result, the az- 
imuthal component of the hydrodynamic interaction vanishes 
to dipolar order. We show that such an azimuthal component 
of the interaction among active rotors arises nevertheless to oc- 
tupolar order and remains finite upon cycle- averaging. While 
the dipolar azimuthal interaction among passive swimmers de- 
cays as 1/r^, the corresponding interaction among active ro- 
tors decays as 1/r^. 

It is important to stress that the following results apply for 
a general active rotor, regardless of its shape : 

(i) The dipolar flow field has no azimuthal component. This 
result is not new and has been known for a long timeB^. It 
is due to the torque-free nature of the swimmers and the fact 
that the rotletlstresslet decomposition coincides with the ra- 
dial/azimuthal one in the plane orthogonal to the axis of rota- 
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tion. 

(ii) The quadmpolar flow vanishes when averaged over a 
period of rotation. This rehes on the two-dimensional nature 
of the problem and the equivalence between the period average 
and an orientation average, which holds if the period of rota- 
tion is much smaller than any other time scale in the problem. 
The vanishing of the cycle-averaged quadrupolar contribution 
then follows immediately from the symmetry properties of the 
multipole expansion of the flow field. 

(iii) The octupolar flow field remains finite upon cycle aver- 
aging and results in an azimuthal interaction force that decays 
as 1/r^. 

Furthermore, these results can be extended to a more gen- 
eral form of the Green's function. We have limited ourselves 
to a discussion of the flow field in planes normal to the rotor's 
axis of rotation. A motivation for this is that rotors are often 
confined to a liquid/air, liquid/liquid or liquid/solid interface. 
On the other hand, such confinement often changes the form of 
the Green's function. Assuming that the medium is isotropic 
in the plane perpendicular to the rotors' axis of rotation, and 
to which they are confined, the most general Green's func- 
tion describing the flow field in this plane can be written as 
^ij{r) =A{r) 5ij-\-B{r) fffj where A and B are two unknown 
function of r = |r|. In addition, we assume that at distances 
large compared to the thickness of the interface ^/j(r) decays 
as 1/r^, i.e. ^ij{r) = {a 5ij -\-b rirj)/r^ where a and b are 
two unknown constants. The three-dimensional Oseen tensor 
used in this work corresponds to ^ = 1 and a = b = 1. The 
Green's function describing interactions in a thin fluid film of 
two-dimensional viscosity T] , surrounded by by a bulk fluid of 
viscosity r]^ corresponds ton= I, a = and b = 2r{/r{e^^ • A 
higher value of n can be found, e.g., if the rotors are close to 
a solid wall with no-slip boundary condition. The 1/r term 
is then canceled by the presence of a virtual source of oppo- 
site s treng th located on the other side of the wall, leading to 
n = 2l ^^l^^ [ The flow field resulting from any of these Green's 
functions is qualitatively the same as the one derived in this 
work. More specifically, the points (i), (ii) and (iii) above re- 
main true regardless of the values of a, b and n, except that the 
octupolar term now decays as l/r"+^ *. 

On the other hand, the exact form of the azimuthal octupo- 
lar interaction between active rotors depends on the nature of 
the rotor. The one-bead swimmer we described in this paper 
is a simple model system in that it has an isotropic drag (via 
the spherical bead) and rotates around itself. A more general 
model would feature an anisotropic drag (e.g. because of a 
non-spherical bead or the presence of more than one bead) 
and/or a significant unperturbed trajectory radius Rq. The az- 
imuthal interaction between two rotors would then depend on 
their relative orientation, the dynamics of which would have to 

* Note that there is a special case, namely na + b ^0, in which there is no 
azimuthal field at all, regardless of the nature of the rotor. 



be worked out. This is beyond the scope of the present paper 
and left for future work. 



Appendix A: Period averaged flow field 

The flow field generated by a rotor depends in general on its 
orientation. As the particle rotates and moves through the 
fluid, the flow field evolves in time as well. This evolution 
is characterized by two well separated time scales: the period 
of rotation Tq of the rotor and the (generally much longer) time 
scale associated with its motion through the fluid due to inter- 
actions with other rotors. On time scales large compared to 
7b, the swimmers' dynamics is governed by an effective flow 
field obtained by averaging the instantaneous flow field over a 
period of rotation. More precisely, the corrections to the av- 
erage flow field u arising from instantaneous fluctuations over 
time scales Tq are of order 5u ^ TquVu ^ u^To/r where Tqu 
is the typical distance travelled by the rotor during the time Tq 
by following the flow field of other rotors located at distance 
r. The leading term of the flow field is at most dipolar, with 
u ^ 1/r^. The correction is then 5u ~ 1/r^ and is negligible 
compared to terms up to octupole order (~ 1/r^) in the aver- 
age flow field. 

Similarly, the fluctuations in the angular velocity are gov- 
erned by the ratio of the curl of the flow field in Eq. ( [4b| ), 
which is at most of order F /r]r' (corresponding to a dipolar 
flow), to the angular velocity (0 ~ Fl/C^^ given in Eq. ( |9b| ). 
As a result, we can assume that the rotation of an individual 
swimmer takes place at a constant speed (up to corrections of 
order 1/r^) and replace the average over a rotation by an aver- 
age over the orientation of the swimmer. 

These observations greatly simplify the dynamics. In par- 
ticular, the quadrupole term in the multipole expansion of the 
flow field, which changes sign when the rotor is rotated by 7t, 
does not contribute to the averaged flow. It should be stressed, 
however, that the angular averages washes out the angular de- 
pendance of the flow field, yielding an isotropic effective flow. 

Finally, a study of the corrections due to these fluctuations 
on the dynamics of passive rotors can be found in Refs. [27ll33l 
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